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c5 '. Abstract 

Within the frame of quantum optics we analyze the properties of spontaneous emission of two- 
level atom in media with indefinite permittivity tensor where the geometry of the dispersion relation 

Q ■ is characterized by an ellipsoid or a hyperboloid (hyperbolic medium). The decay rate is explicitly 
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pi |. given with the orientation of the dipole transition matrix element taken into account. It indicates 
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^ ' that for the ellipsoid case the intensity of the photons coupled into different modes can be tuned 
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' S ' ^y changing the direction of the matrix element and for the hyperboloid case it is found that spon- 
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background. Moreover, spontaneous emission exhibit the strong directivity and get the maximum 
in the asymptote direction. 
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I. INTRODUCTION 



Metamaterials is composed of the periodic dielectric element arrays where metal mico- 
structure e.g.the split-ring resonators(SRRs) is included in each cell and used as a effective 
continuous medium in a narrow microwave frequency region [l|-y] or in the near- visible light 
regionjSi]. The materials have the left-handed property for its negative permittivity and 
permeability simultaneously and cause to the refocusing and phase compensation [6|, which 
provide a new manipulating space for the design of quantum optical devices. For instance, 
the suppression of spontaneous emission and superradiance over macroscopic distances in 
the left-handed materials (LHM)[7|, the quantum interference enhancement between two 
spontaneous emission transitions with the LHM[8j], long-lived entanglement between two 
distant atoms via the LHM[9i]. On the other hand the effective permittivity or permeability 
of the metamaterials only holds for the specific Ipolarization and it mean that metamaterials 
is anisotropic as well as left-handed 10l-ll2|. The anisotropy of metamaterials has also been 
found some new optical properties for example, the incident electric field E can couple to the 
magnetic resonance of the SRRs when the electromagnetic waves propagate perpendicular 
to the SRR plane 13|] and it can be utilized to excite electron-spin resonance 14l |. 



Since the elements of permittivity(permeabi 



ity) tensor of the metamaterials can be neg- 



ative or positive in different frequency range jlO| we apply the term indefinite to anisotropic 
media in which not all of the principal components of the "?" tensor has same sign. The ge- 
ometry of the dispersion relation in the indefinite media can be characterized by an ellipsoid 
or a hyperboloid as shown in the inset to FiglH The medium with hyperboloid geometry of 
the dispersion relation is also called "hyperbolic medium" (HM). The HM has been recently 



found to be having many novel properties such as the superlens effect |15Nl7|. slow-light 



effect |18l| and the "big flash" of the photons in the HM by an optical metric signature phase 
transition [19(]. Moreover, some interesting quantum optical properties(QOPs) of the HM 



has also been found experimentally 201] and theoretically 2ll-l23|. for instance, controlling 
spontaneous emission with the HM[2C!|], broadband Purcell effect 211]. the dipole radiation 



and its enhancement near the surface of the HM|22|. l23l]. There are also some studies on the 
QOPs of the indefinite media where the emphasis is put on the left-handed property[24] and 
the singularity of the density of states [25]. 



The above theoretical analysis on the QOPs of the indefinite media are within frame- 




FIG. 1: The schematic diagram of the three vectors, inset: The geometry of the dispersion 

relations of the indefinite medium 



work of the macroscopic electromagnetic wave theory and the atom is approximated as a 
dipolejSl, I2ll-l23( where the radiated power is exacted from the Green's function of the sys- 
tem. However, the rate of spontaneous emission for a two-level atom in the indefinite media 
has not been investigated within the framework of the quantum optics. In this letter we will 
explicitly give the expression of the decay rate for a two-level atom in the indefinite media 
under the Weisskopf-Wigner approximation. It indicates that for the ellipsoid case the in- 
tensity of the photons coupled into different modes can be tuned by changing the direction 
of the matrix element. For the hyperboloid case it is found that spontaneous emission in 
the HM can be dramatically enhanced in comparison with the dielectric background, mean- 
while, the spontaneous emission exhibit the strong directivity and get the maximum in the 
asymptote direction. 



II. MODEL AND FORMULA 



To deal with the our problem some theories of the quantum optics of dielectric media 
is needed. There are many schemes of the electromagnetic quantization for the different 
types of media. One of them is that the quantities characterizing the dielectric such as 
polarization field, magnetization field or some quantities of their combination are involved 
in the quantization procedure and the interaction between their quanta and the photon 



is taken into account 26N28l|. Another scheme is that the permittivity and permeabihty 
tensor are regarded as the parameters in the field equation with different gauges (gauge 



conditions) chosen for the different media 



29 



3l| . In our case the medium is considered 



as anisotropic, homogenous, non-dispersive and lossless for simplicity. The quantization 
scheme in Ref. 3l| is used where the medium is characterized by the constitutive equations 
D(r) = W\r) ■ E(r) + "e^^^r) ■ B(r); H(r) = "fl^^^r) ■ E(r) + p'^\r) ■ B(r). For our case 
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where e^ > and e^ > for the uniaxial anisotropic medium or e^ > and e^ < 
for the HM. According to the Maxwell's equations the dispersion relation for the medium 
characterized by Eq.([l]) in principal axis coordinate system is expressed as: 






(2) 



The geometry of the Eq.(l2]) in the K space represent a hyperboloid or an ellipsoid that 
depends on the sign of et where e/, > is assumed. Here we are interested in the decay rate for 
a two-level atom in the indefinite media. As the excitation of the quantum electromagnetic 
field the photon is emitted form an atom and subsequently coupled into the field modes 
permitted by the dielectric environment. In dielectric system the photon propagates in the 
fashion of the classical field modes of the system and inversely interact with the atom. In 
the spirit of Einstein's original model the total energy of the electromagnetic field mode 
in the dielectric should be ^fcJ32[, which can be used to determine the field amphtude in 
the factor gk- According to Ref. 3l| the eigenvector field in the dielectric should satisfy the 
following gauge condition and the eigenequations: 
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t(r) ■ Fk(r) 
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The quantized electric field can be expressed as E{r, t) = ^^ ekFk{r)e *'^'=*dfc + H.c. and the 



total Hamiltonian of the photon and the atom under the rotating-wave approximation is 

k 

+ ^ ^ gk{cr+ak + <3-_4) (5) 

k 

where dz = \a){a\ — \b){b\, &+ = \a){b\, (5"_ = \b){a\, \a), \b) is the excited and ground states 
of atom with eigenvalues Ea, Eb and 

9k = ^ (6) 



^ab is the matrix element of atom's dipole between states ^ab = —{a\eR\b) = —{ 



eR\a). 



For simplicity the Hamiltonian Eq.(l5]) can be expressed as in the interaction [33|] picture 

k 

where Ea — E^ = fiw^. The state vector of the composite system of the photon and atom 
including the vacuum state |0) is IV'(^)) = Ca|a,0) + 'Yl,k^b,k\bAk)- With the Weisskopf- 
Wigner approximation the decay rate of the atom in the anisotropic medium can be expressed 
as the following integration 
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where the dispersion relation of the modes is written as, u = f{k, 6, 0) and A; is A; = 
f^^{uj,9,(l)). Here the eigen electric field is assumed as Ek{r) = Fk{r); H^ = ^^V x E^. 
That a single photon is coupled into the eigen-mode require the total energy of the mode 
U = ^ / [^k -^ ■ Ek + fio\Hk\'^] d^x = hu, which determine the amplitude Ffc(O) under the 
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box normalization condition. 



A. Spontaneous decay rate in the medium w^ith dispersion geometry of ellipsoid 

In order to check the validity of the formulism, the case e^ > is firstly considered 
where the medium become uniaxial anisotropic with an ellipsoidal dispersive geometry. In 
the medium there are usually two types of eigen-modes named as extraordinary wave and 
ordinary wave which accommodate the emitted photons form the atom. To this end we have 
to explore the amplitudes and the energy of the two modes in detail. The eigen electric field 



is supposed as the plane wave Ek{r) = E^qc^^'^ e^'^^^^^ and Hk{r) = Hko^^'^'^ e^^^^^^ where the 
complex vector amplitudes are to be determined. Substitute the expression into the Maxwell 
Equations and after some algebra we get the two sets of dispersion and polarization relations 
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In the spherical coordinate system the dispersion relations can be uniformlly written as: 

■j= for the transversal mode EqQ, and f{6) 



u{k) = f{9)ck, where f{9) 
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for the lognitudinal mode Eq( !TO|) . For the transverse mode the total energy in the volume 
V and the relative amplitude are 
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For the lognitudinal mode the corresponding quantities : 
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In the principal axis spherical coordinate system there are three vectors to be identified. 
They are the vectorial transition matrix element ^ab = i^ab, ^'o, ^o)) the electric field vector 
Eko = [Eko, Oi, 0i) and the wave vector k = [k, 6, 0) which is shown in Fig. [H The factor in 
Eq.dH]) is explicitly given 

, |2 ^IbElpCOS^Oo,! ..o^ 

l^^l = ^ '-1'^'' 

where ^o,i is the angle between !^ab and Eko- According to the geometrically relations of the 
vectors ^ab and Eko and the gauge condition k -^ ■ Eko = , we get the equations cos 6^0,1 = 
sin 9o sin 6*1 cos(0o — 0i) + cos Oq cos 61 and ei cos 6*1 cos 9 + et sin 9i sin 9 cos (0i — 0) = 0. In 
term of Eq.(l9l).Eq. (iT0l) . It is noted that 0i = + f for the transversal mode and 0i = for 



the lognitudial mode. With these conditions we get the factor c, = cos^ 6'o,i for the different 
modes 

cJ^ = sin2^osin2(0o-0) (14) 

^ _ [et cos 6q - ei sin Op cot 9 cos {(f)o - (f))f 
""« ~ eicot^^ + e? ^^^^ 



After the implementation of the Eq.Q we get the decay rate for the two modes 
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When eL = et > 0, Eq. fll6cp reduce as 3°^|l°g ;-f = r(et)) which give the decay rate of the 
atom in the homogeneous isotropic medium with permittivity et- The decay rate relative to 
the vacuum T/Tq = F in the case et = 2.5, ej^ = 3.5 is shown in Fig|2]for the different modes. 
Form the Figl2] it is found that when the matrix element vector ^ab is parallel to the k^ 
axis(^o = 0), the atom can get the maximal coupling with the lognitudial mode of the system 
where Ez y^ and much more photons emitted form the atom is coupled into the mode. 
Meanwhile, the transverse mode get the no coupling with S'ab for Ez = and Tt = 0. In the 
case, F = F^, = 1.58 = y/e[, which the anisotropic medium behave as the isotropic medium 
with permittivity et for the atom's spontaneous emission. With the increase of ^o, Tt also 
increase due to the enhancement of the coupling with the transverse mode and Tl decrease 
due to the reduction of the coupling with the lognitudial mode. When 6*0 = f where ^ab 
is perpendicular to the kz axis Ft get the maximum and F^ get the minimum. According 
to Eq.f llGcp F get the maximum(et < e^,) or the minimum(ei > e^): F^ = ^^^°3^^^^ ^ at 
6*0 = f. These results can be used to control the intensity of the different modes from the 
spontaneous emission by tuning ^o- 

B. Spontaneous decay rate in the medium with dispersion geometry of hyper- 
boloid 

For the case et < 0,6^ > 0, Eq.([l]) indicates the hyperboloid geometry of the dispersion 
relation. Since et < 0, the branch k = ^a/c* can only exist with evanescent wave. This 
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FIG. 2: The decay rate of atom in medium with e^ = 2.5, ez, = 3.5 relative to the case for 

the vacuum F = F/Fq. Doted hue: the decay rate for the transverse mode. Dashed hne: 

the decay rate for the lognitudinal mode. Sohd hne: the total decay rate 



mode can not carry the energy away from the atom for the spontaneous emission and the 
contribution to the decay rate is ignored for simplicity. For the branch ^^ — - + -^ = (7) 
the decay rate F is proportional to the following expression: 



Y = fiC dx^^^^' ^' " ^ ^"^^ " ^^ '* '°'' ^° 



where x = cos 9, e^ 
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It is noted that under the Weisskopf-Wigner approximation the integration Eq.(|8]) is 
actually calculated on the equal- frequency surface. For the ellipsoid Eq. ([8]) 9 integrates over 
the range [0, vr], while for the hyperboloid 9 integrates over the range [9a, tt — 9a], as shown 
in FigJ3] where 9a is the polar angle of the asymptote. In term of Eq.f lT7|) there are two poles 
±V^ the positions of which on the axis depend on the parameters et, ei- When e^ > e^ > 0, 
e^ > 1 and ei > tt > 0, e^ < the poles ±A/e^ are out of the range [—1, 1] or on the imaginary 
axis. It enable Eq.f lT7|) to be calculated and the result is given in Eq.f lT6|) . However, for 
the hyperboloid case where et < 0, e/, > 0, < e^j < 1 the poles lie in the range [—1,1] 
the integration diverges. This divergence is the manifestation of the change of the topology 
from the ellipsoid to the hyperboloid and it also indicate that hyperboloid of the dispersion 
relation is only the perfect effective medium approximation of some composite materials e.g. 
photonics crystal, metamaterials under the long wave-length limit. Therefor, some cutoff 
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FIG. 3: The schematic diagram of integration region for the hyperbohc medium where 

et < 0,eL > 



methods have to be introduced for the calculation of Eq.( fT7l) in the case et < 0, e^ > 0. To 
this end it is defined that Xa = cos 6a = y^ and cos^c = Xc = axa, where ±Xc is the new 
integration limits for Eq.( TT71) with condition Xc < Xa- The integration range of the variable 
6 for the cutoff is marked in the shadowed region in FigJS] When the coefficient a = 1, the 
integration limit approach to the two polar poles ixa- Under the cutoff approximation we 
get the decay rate T as follows: 



T. T- "^ [^i - ^L cos^ 6^0 i-Aeiet + e\ + 6e?)] - Qaeiet cos^ Oq (ez, - e*) 

^H = ^o — -vi:r. ~ — ; Uoj 



Aia^-l 
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eL - Q 
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where Fn 



^ab^l 



In FigJH we give the Th = F/Fq vs the parameter a for the different 6*0 with e^ = 
—2.5,eL = 3.5. It is obviously found that Th increases on the whole when a approach 1. It 
is noted that when a close 1 the Th increases more sharply. It is because that more large 
a is, more modes with high k vectors are involved in the spontaneous emission. In HM the 
modes with high k points usually have large density of states and more photons emitted 
from the atom can be accommodated. Moreover, since the photons is mainly coupled into 
the lognitudial mode, Th get more large value when the matrix element vector !^ab get the 
small angle Oq with z axis which enhance the coupling factor g^. 
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FIG. 4: The relative decay rate of atom in medium with ej = —2.5, ei = 3.5 relative to the 
case for the vacuum according to parameter a for different angles ^o 
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FIG. 5: The relative decay rate of atom in medium with e^ = —2.5, cl = 3.5 to the case for 
the background e^ = 3.5 according to parameter a for different angles 6o 

To explore the effect of the enhancement of the spontaneous emission we compare the 
decay rate T^ with the decay rate F^^, the decay rate in the background medium with 
permittivity e^. The relative decay rate Th = Tfj/^eL corresponding to parameter a is given 
in FigEl It is easily found that T^ can be lager than one only if a is sufficiently large in spite 
of the difference of ^o, which mean that the more like the strict hyperbolic medium the real 
materials behave, the more the spontaneous emission is enhanced relative to the background 
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FIG. 6: The amplitude of the integrand 7(6*) corresponding to 9 with parameters 

e^ = -2.5 + 0.5^,6^ = 3.5 



medium. According to Eq. (IT8ll there are IrniT h = c>o, which imphcate that the perfect 
hyperbohc medium is the an ideahzation model of some real composite materials. Besides 
the enhancement of the spontaneous emission the directivity is also worthy of being noted. 
To this end the Eq.( !T7|) can be rewritten as F = (3 f^ ■y{9)d9. Considering the divergence, 
a small imaginary part is added to e^. Figj6] gives the amplitude of the integrand 7(6') with 
et = —2.5 + 0.5i, e^ = 3.5. There are two obvious peaks in the positions corresponding to 
the directions of the asymptotes. This strong directivity of the spontaneous emission has 
been noted experimentally and used to design the single gun | la. l34 |. 



III. DISCUSSION 



We explore the spontaneous emission of the two- level atom in the homogeneous anisotropic 
medium where the dispersion geometry exhibit as an ellipsoid or a hyperboloid and the cor- 
responding decay rate F in detail under the Weisskopf-Wigner approximation. For the 
ellipsoid case, there are two kinds of modes that contribute the decay rate F and to some de- 
gree the medium with ellipsoid dispersion provides two type of mode space to accommodate 
the emitted photons. Moreover, the polar angle 9o also be used to change the intensity of 
the photons coupled into the different modes. When ei = e^ the obtained formula Eq. fll6cp 
reduce to the case of the isotropic medium, which justify the validity of our model. 
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FIG. 7: Left: Schematic view of the "layered" metamaterials. Right: the corresponding 

equvi-frequency contour 



When the above model is applied to the hyperboloid case the divergence is encountered 
and the cutoff for the integration variable 6 (equivalently, the wave vector K) is introduced 
for the approximation of the real materials. Though the transverse radiating wave mode 
degenerates into the evanescent wave mode, the enhanced spontaneous emission relative 
to the background medium e/, can be obtained only if the dispersion geometry of the real 
materials is sufficiently close to the strict hyperboloid. It mean that more high k modes which 
corresponds lager parameter a get involved in the coupling of the photons. In real world, the 
hyperbolic medium is used as the perfect model on the some materials with periodic micro- 
structure where band structure exhibit the approximative hyperbolic geometry in the small 
range of fc, e.g. the hyperbolic metamaterials. The simplest of them is the one dimensional 
photonic crystal including the metal layer as shown in Fig. [7]together with the equifrequency 
contour. Generally, in the long wavelength limit the system can be approximated as the 
hyperbolic medium within the effective medium theory. It is expected that when the more 
small the dimension of the lattice a is, the more close to the perfect hyperbolic medium the 
system is and more enhancement of the spontaneous emission can be achieved. 

The physical mechanism of the enhancement on the spontaneous emission in HM can 
be understood from the perspective of density of states. According to the Fermi's golden 
rule the one-to-many transition probability per unit of time depends not only on the matrix 
element but the density of final states p{uj) as well. When the geometry of physical dispersion 
relation of the medium changes form the ellipsoid to the hyperboloid, due to the change of 
the topological property the density of states diverges in the lossless continuous hyperbolic 



medium limit: p{uj) 
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where Krvt is the momentum cutoff 
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is defined by either metamaterial structure scale or by losses. It is the occurrence of the 
large p that enable the transition probability to be increased dramatically. Even the loss and 
the dispersion of the e{(jj) in the materials is taken into account the spontaneous emission is 
expected to be largely enhanced. 
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